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ABSTRACT 

Commissioning observations with the Apache Point Observatory Galactic Evolution Experiment 
(APOGEE), part of the Sloan Digital Sky Survey III, have produced radial velocities (RVs) for ~4700 
K/M-giant stars in the Milky Way bulge. These high-resolution (R ~ 22,500), high-S/N (>100 per 
resolution element), near-infrared (1.51-1.70 /im; NIR) spectra provide accurate RVs (ev~0.2 km s _1 ) 
for the sample of stars in 18 Galactic bulge fields spanning — 1°< I < 20°, \b\ < 20°, and 5 > —32°. 
This represents the largest NIR high-resolution spectroscopic sample of giant stars ever assembled in 
this region of the Galaxy. A cold (av~30 km s _1 ), high- velocity peak (Vgsr~+200 km s _1 ) is found 
to comprise a significant fraction (~10%) of stars in many of these fields. These high RVs have not 
been detected in previous MW surveys and are not expected for a simple, circularly rotating disk. 
Preliminary distance estimates rule out an origin from the background Sagittarius tidal stream or a 
new stream in the MW disk. Comparison to various Galactic models suggests that these high RVs are 
best explained by stars in orbits of the Galactic bar potential, although some observational features 
remain unexplained. 

Subject headings: Galaxy: bulge — Galaxy: kinematics and dynamics — Galaxy: structure — surveys 



1. INTRODUCTION 

Multiple Af-body models have demonstrated that ro- 
tating stellar disks form bars rela tively quickly under 
a wide var i ety of conditions (e.g., Bureau fc Athanas-| 
soula 2005; Deba ttista et al.||2006 |Martmez-Valpuesta| 
et al. 120061 , and observations nave identified a high inci- 
dence ol Targe centra l bars in extragalactic sys tems (e.g., 
Liitticke et aLl|2000j |Chung fc Bureau||2004| ). Even in 
edge-on spiral galaxies, where the stellar bar can be dif- 
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ficult to isolate photometrically from the surrounding 
bulge and/or disk, the bar can be identified via kine- 
matical signatures: double-peaked rotation curves, radial 
velocity (RV) dispersion profiles with particular combi- 
nations of plateaus and peaks/minima, and the relation 
between the mean and higher moments of t he RV distri- 
bution (e.g., Bureau & Athanassoula 2005). 



In recent decades, good evidence tor a Milky Way 
(MW) bar has arisen from infrar ed (IR) surface bright- 
ness maps and star counts (e.g., |Blitz fc Spergel|[ l991; 
Weinberg 1992; Dwe k et al.||1995[ |H ammers ley et al . 
2000 ; ICole fc Weinberg! |2QQ2nBen.iamin et al||2005 



Lopez-Corred oira et al. 2007 MRobin et al.|2012p 7 carbon 
stars ( |Cole fc We inberg 2002 ), the solar neighbo rhood 
velocity held (e.g.,[Dehnen 2 000] |Minchev et aLp007| ), 
microlensing (e.g., |Gyuk|P^99[ jKattenbury et al. [2UU7 ), 
and inner Galaxy gas kinematics (e.g., [Liszt fc Burton 
1980 |Binney et aLl[l99T| [Fux1[l999i | Weiner & Sellwood 
1999). However, the high extinction towards the inner 
Galaxy midplane (Ay «5-10) has long precluded an ex- 
tensive spectroscopic survey of the type required to mea- 
sure the bulge/bar stellar kinematics on par with studies 
of extragalatic bars in galaxi es viewed edge-on . 

Mor e recently, the BRAVA ([Rich et al |2007b[ ) and AR- 
GOS (Ness & Freeman 2012) surveys have begun per- 
forming kinematical measurements in the bulge periph- 
ery, mapping the bulk cylindrical (i.e., bar-like) rotation 
of the bulge. In this Letter, we report the discovery, 
within data from the new APOGEE survey, of a "cold" 
kinematical feature in the inner disk, which has not been 
seen in previous surveys and is, most likely, a signature 
of particular stellar orbits in the MW bar. In ^2j we 



2 



NIDEVER ET AL. 




Fig. 1. — (a) Coverage of APOGEE commissioning fields (white circles) in the Galactic bulge and inner disk regions. The eight fields with 
distinct high- velocity peaks are marked with black dots. The background image is the 2MASS K s band (Skrut skie et al.||2006| > showing 
the boxy bulge. The Sagittarius dwarf spheroidal is indicated by an ellipse, (b) Reddening-corrected Hess diagram ot all observed stars 
(grayscale) with the high-velocity stars with [J — K s ]o > 0.5 shown as red dots. The high- velocity thresholds in Table [l] which represent 
the RV of the "trough" betw een the main peak and the high-velocity peak, were used to select the high-velocity stars from each field. 
Fiducial Girardi et al. (2002} isochrones (2 Gyr, [Fe/H]=0.0) are shown at distances of 5/10/30 kpc (blue/green/red). 



briefly describe the observations and data reduction, and 
we present in ^3] the RV distributions of the APOGEE 
fields, along with the RV predictions of multiple Galac- 
tic kinematical models. In Q we use these comparisons 
to interpret the nature of the APOGEE observations, 
and we evaluate alternative explanations for the data and 
model discrepancies. 

2. OBSERVATIONS AND DATA REDUCTION 



The APOGEE project ( |Eisenstein et al.||2011[ Majew- 
ski et al. in p reparation) uses a custom-built, cryogenic 
spectrograph ( | Wilson et al.|201Q[ ) recording 300 simulta- 
neous near-IR spectr a"? 1.51-1.7U£mi ) fed from the Sloan 
2.5-m telescope ( |Gunn et al. 2006| ). The final config- 
uration for the spectrograph is described in Wilson et 
al. (in preparation); however, the data described here 
were taken before the instrument was in its fully com- 
missioned state and in optimal focus (achieved Septem- 
ber 2011). The net effect was a somewhat blurred line- 
spread-function, yielding a degradation of the resolution 
primarily in the red detector to R ~16,000 (1.65-1.70/im) 
compared to R ~ 23,000 in the other two detectors 
(1.51-1.58/im and 1.59-1.64/im). 

The APOGEE fibers are plugged into standard Sloan 
2.5-m plugplates and observed similarly to those of the 
optical Sloan spectrographs, with the following varia- 
tions: (1) 35 fibers in each plugplate configuration collect 
sky spectra, (2) another 35 fibers are placed on bright, 
hot stars to gauge telluric (CO2, H2O, and CH4) absorp- 
tion, and (3) because the large zenith distances required 
to observe the bulge from APO cause strong differen- 
tial refraction, the field diameter was restricted to 1° for 
/<12° fields (2° for the rest). 

We observed 18 bulge fields in June- July 2011 with 



-1°< I < 20°, \b\ < 20° and S > -32° (see Fig. Fiji). 
Targ ets were selected from the 2MASS Point Source Cat- 
alog ( |Skrutskie et al.|2 006) requiring (J — K s )o > 0.5 (al- 
though some plugplates had no color selection for testing 
purposes) and i^<11.0 (see Zasowski et al., in prepara- 
tion, for details on APOGEE targeting strategy) pj The 
dereddened Hess diagram for all observed stars is shown 
in Figure [TJ). Each field was observed for a total integra- 
tion of roughly one hour, for a S/N>100 per APOGEE 
pixel for each star. 

An automated data reduction pipeline written specifi- 
cally for APOGEE data (Nidever et al., in preparation) 
was used for (1) collapsing the sample-up-the-ramp dat- 
acubes to 2D images with rejection of cosmic rays, (2) 
extraction of ID spectra from the 2D images, (3) wave- 
length calibration using ThArNe/UrNe arc lamp expo- 
sures and airglow lines, (4) correction of telluric absorp- 
tion as monitored by the 35 hot star standards, (5) sub- 
traction of airglow lines and sky continuum, and (6) 
derivation of RVs by cross-correlation against a grid of 
template spectra spanning the expected range of temper- 
atures, gravities, and metallicities of target stars. The 
median RV scatter for stars observed multiple times is 
0.22 km s _1 (with a mode of ~0.09 km s _1 ), while tests 
of derived RVs for observations of 53 stars in three glob- 
ular clusters (M3, M13, M15) with accurate literature 
RVs suggest that the zeropoint of the APOGEE pipeline 
RVs is accurate to -0.26±0.22 km s" 1 . 



15 Th e 2MASS photometr y was deredd ened using the R JCE 
method (Majewski et al. 2011 ), with WISE ( [Wright et al.|2010| > and 
GLIMPSE (|Churchwell et al. 2009) providing the required mid-IR 
photometry. 
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Fig. 2. — Velocity (Vqsr) histograms of eight APOGEE bulge fields showing a dual-peak structure, using 20 km s _1 binning. A two- 
Gaussian model (see text) is fitted to the data to determine the central velocity of both peaks (blue/red lines). The solid black line is the 
sum of the two components. 



3. 



RADIAL VELOCITY DISTRIBUTIONS: DATA AND 
MODELS 



The bulge RV distributions exhibit a significant (and 
unexpected) number of high- velocity stars, which appear 
as prominent high-velocity peaks in eight fields and as 
high-velocity "shoulders" to the main peak in many of 
the other fields. We focus on the RV distributions of the 
eight APOGEE bulge fields with distinct high-velocity 
peaks, shown in Figure|2]with a (J—K s )o > 0.5 selection. 
The stars in the high-KY peak account for ~10% of the 
stars in a field. To derive the mean velocity of the high- 
velocity component, the histograms are fitted with a two- 
Gaussian model - a hot component for the bulge and a 
colder component for the high-velocity starsj^] 

These high-velocity peaks have not been previously 
detected in the MW bulge, including in the ext ensive 
BRAVA data (-10°< I < +10°, b = - 



-6°, -8° 



Kun- 



der et al. 2012). Several narrow peaks in t he RV his 



togram w ere detected in the BRAVA data (Howard et 
|al.||2008) but were eventually determined to be due to 
stochastic effects. Our high- velocity peaks are unlikely to 
be due to stochastic effects because: (1) our sample per 
field is large (twice that of BRAVA), (2) the high- velocity 
peaks are seen in at least eight APOGEE fields, (3) the 
peaks appear at similar RVs following a general trend 

16 The models are not intended to be perfect "physical" rep- 
resentations and are poor fits for several fields (e.g., [4.0,0.0] and 
[5.7,-2.0]), but they serve well to estimate the peak center. 
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and (4) the pattern is "mir- 
rored" across the miclplane (e.g., strong similarities be- 
tween Z, 6=10°, ±2°). We conclude that the high- velocity 
peaks are due to real motions of stars in the Galaxy. 

Table [I] gives the field name, longitude, latitude, di- 
ameter, mean Vqsr and velocity dispersion of all stars, 
number of stars, mean Vqsr of the high-velocity stars, 
the high- velocity threshold (essentially the velocity of the 
"trough" between the main and high- velocity peaks) , and 
the number of high- velocity stars in the eight APOGEE 
bulge fields that exhibit the most distinct high-velocity 
peaks. To ascertain the nature of the high- velocity stars, 
we compare the APOGEE RV histograms to several 
models. 



3.1. Besangon Galaxy Model 
The Besancon Galaxy Model dRobin et al.| 2003 



2012 



hereafter BGM) simulates the stellar content in any given 
line of sight and computes the photometry, kinematics, 
and metallicity of each simulated star. For each popula- 
tion (thin disk, thick disk, halo, bar, and bulge), a star 
formation rate history, age, and initial mass function are 
assumed, allowing the model to generate the distribu- 
tion function of absolute magnitude, effective tempera- 
ture, and age of the stars. Density functions are assumed 
for each population and tested against observations using 
photometric sta r counts. The model also includes a 3D 
extinction map ( |Marshall et al.||2006| ) and simulated un- 




Fig. 3. — (a) Position-velocity diagram of the two-component Gaussian fits to the APOGEE bulge RV histograms, with main peaks in 
blue and high-velocity peaks in red. Gaussian FWHM values are indicated by the vertical lines. A general trend is seen in the high-velocity 
component, with mean velocity first rising with longitude and then slowly falling at longitudes above 1^6—7°. The mean velocities of the 
main components fall close to the theoretical disk trend calculated from a simple rotation curve (Vqsr— 220 sin/; black), (b-e) Spatial 
distribution of particles in the Kazantzidis 200 pc scale-height model (upper panels, b+c) and the MVG model (lower panels, d+e). All 
particles are shown in the left panels (b+d) while the right panels (c+e) show the distribution of particles with Vqsr >150 km s _1 which 
live on the leading edge of the bar. The sun's position is at (0,-8.5) which is off the lower edge of the panels. The black contours outline 
the general shape of the model MW bar. 



certainties on the observational parameters of each star. 

The BGM bar a nd bulge populations are described in 
Robin et al.| ( |2Q12[ ) . The bar is the more massive compo- 
nent and dominates the stellar content at low latitudes, 
while the bulge is longer and thicker and gives a contribu- 
tion at intermediate latitudes where the bar becomes less 
prom inent. The bar kinematics are taken from the Fux 
( |1999| ) dynamical model, and the bulge kinematics are es- 
tablished to reproduce the BRAVA survey data (Robin 
et al., in preparation). This multi-population model 
also explains well the presence of double red clumps at 



med ium l a titudes (|Nataf et al. | [201Q[ |McWilliam fc Zoc- 
cali 2010j Saito et al. ' 2Ullj ) that are created by stars 



trap ped in vertical reson ances associated with the bar 
(e.g., Combes et al.|1990 ), lifting them to high altitudes. 



3.2. Kazantzidis N-body Model 

We also co mpare to the A^-body model of Kazantzidis 
et al. (2008), which represents a MW-sized galaxy with 



initial disk scale-height of 200 pc. The model was sub- 
jected to a cosmologically- motivated satellite accretion 
history and formed a bar as a consequence of the accre- 
tion events. The properties of the bar are similar to those 
of the MW bar (in overall orientation and size). 

3.3. Martinez- Valpuesta & Gerhard Bar Model 
The mod el developed in [Martinez - Valpuesta fc Ger- 



hard (2011, hereafter MVG) is an Af-body galaxy based 



on secular evolution of disk galaxies through angular mo- 
mentum transfer between disk and halo. Once the bar 
appears, it grows stronger and buckles, forming a boxy 
bulge, such as the one seen in the MW. After this buck- 



ling the bar weakens and resumes growth by giving an- 
gular momentum to the live dark matter halo. At this 
point the bar has a short, thick component and longer, 
thinner component. This model has already been quanti- 
tatively comp ared with the innermost region of the MW 



(-8°< I < 8°;|Gerhard fc Martinez- Valpuesta] 



also in the outer bar ( jMart inez- Valpuesta 
2011) with good agreement" 



and 



Gerhard 



3.4. Comparison to Models and Other Data 

Figure H] compares the observed APOGEE RV his- 
tograms (black solid lines and shaded region) to the 
above models. 

The blue line shows one "observation" of the BGM, se- 
lected using the simulated photometric information and 
APOGEE targeting algorithms with model photomet- 
ric errors and no reddening. The BGM reproduces the 
main features of most fields quite well, but not the high- 
velocity peak. Multiple simulated observations of the 
BGM were performed and found to be very similar, in- 
dicating that stochastic effects were small. Only one ob- 
servation is shown for clarity. 

The RV histograms of the Kazantzidis model are shown 
in green. This model cannot be mock "observed" like 
the BGM because synthetic photometry is not available. 
The RV distributions are a fairly good representation of 
the APOGEE histograms, but sometimes show a poorer 
match to the main RV peak, compared to the BGM. This 
model does show significantly more stars at high velocity 
and roughly match the fraction and distribution trend 
(decreasing number to even higher velocities) observed 
in the APOGEE data. However, it does not reproduce 
the "trough" between the main and high- velocity peaks 



APOGEE Bar Kinematics 



5 



o 



O 



O 
OS 



o 



o 

OS 




-200 200 -200 200 

V GSR (km/s) V GSR (km/s) 



Fig. 4. — Velocity (Vqsr) histograms (fractions) of the APOGEE fiel ds fblack solid line and shaded region), compared to a variety of 
models and other data. The p opul ation synthesis "Besancon" model of Robin et al. (2003 2012), which incorporates multiple kinematical 
prescriptions as desc ribed in £]3.1| (BGM, blue). Th e A/-bodv model oifR azantzidis et al. (2UU8) with initi al dis k scale-height of 200 pc 
as described in ^3.2| (green). The TV-body model of Martinez- Valpuesta & Gerhard (2U11) as described in ^3.3| (MVG, red). Integrated 
HI data within the field's FOV from the LAB survey ([^Liberia et al. 2UUo"] dotte d black line scaled up b y a tact or of five). Predicted 
distribution of stars from the Sagittarius dSph streams, according to the model of [Law & Majewski (2010 purple line scaled down by a 
factor of three). The Sagittarius model has been limited to points that could possibly tall withm the APUGEE sample (i.e., we assumed 
all points had tip- RGB magnitudes, calculated their apparent magnitude, and removed all points too faint to be included in the sample). 

seen in the APOGEE data. 

The MVG model (red), as with the Kazantzidis mod- 
els, reproduces the high- velocity stars much better than 
the BGM, although often "over-predicting" them. As in 
the Kazantzidis model, the MVG high- velocity stars are 
in the bar. This model also shows the interpeak "trough" 
in some of the fields (e.g., [14.0,-2.0], [6.0,0.0], [4.0,0.0]). 



The Lei den- Argentine-Bonn (LAB) H I data ( |Kalberla 
|et al.|[2005| , shown as a dotted black line, display high- 
velocity peaks in some of the fields, but often at an RV 
different from the APOGEE peaks. Gas kinematics are 
complicated in the inner Galaxy but are d ominated by 
forces influenced by the bar (e.g., Fux|1999j ). There does 
not appear to be a strong gaseous component of this new 



high- velocity stellar feature. 

Finally, because the Sagittarius dwarf spheroidal 
galaxy (Sgr) and part of its tidal tails lie behind the 
bulge in this genera l direction, we test t he A POGEE RVs 
against the recent Law & Majewski (2010) Sgr model 
(purple). It predicts very tew Sgr particles in our fields 
and at RVs quite different from the APOGEE high- 
velocity peaks. 

4. DISCUSSION 

We detect high- velocity peaks in our APOGEE bulge 
RV distributions previously undetected in the data of 
BRAVA or any other survey; this might be due to the 
different regions sampled by the surveys (i.e., APOGEE 
probing lower \b\). We compare our stellar data to var- 
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TABLE 1 
APOGEE Bulge Fields 



ield Name 


^cen 


been 


Diameter 


(Vqsr) 




Nstars 3, 


(VbsR>high 


Hhresh b 




(cleg) 


(cleg) 


(cleg) 


(km s- 1 ) 


(km s- 1 ) 




(km s- 1 ) 


(km s- 1 ) 


004+00 


4.0 


0.0 


1.0 


39.4 


100.9 


376 


203.6 


141 


004-04 


4.3 


-4.3 


1.0 


47.1 


89.8 


246 


214.5 


156 


006-02 


5.7 


-2.0 


1.0 


39.5 


88.4 


245 


200.0 


146 


006+00 


6.0 


0.0 


1.0 


50.2 


96.1 


364 


234.3 


177 


010-02 


10.0 


-2.0 


1.0 


72.0 


80.6 


249 


221.7 


170 


010+02 


10.0 


2.0 


1.0 


69.5 


76.2 


247 


233.0 


180 


014-02 


14.0 


-2.0 


2.0 


93.7 


77.3 


247 


206.8 


176 


014+02 


14.0 


2.0 


2.0 


85.2 


85.2 


251 


195.1 


171 



Nhigh (%) 



56 (14.9) 

31 (12.6) 
29 (11.8) 
39 (10.7) 
28 (11.2) 
22 (8.9) 

32 (12.9) 
36 (14.3) 



a Some fields have more than 300 stars due to multiple (distinct) plates being observed. 
b The velocity threshold for the high-velocity stars. 



ious models and H I data to discern the nature of this 
newly found "stellar population" of the bulge. 

One potential explanation of this RV feature is that 
these stars are part of the tidal tail of Sgr, which lies 
close to the APOGEE bulge fields and has a m ean ve- 
locity of km s _1 at these longitudes (|Law 



Majewski|2Q10[ ) . However, initial distance estimates (us- 



ing preliminary APOGEE stellar parameters and Girardi 
et al. 2002 isochrones) of -5-10 kpc for the APOGEE 
stars (also see isochrones in Fig. [TJd) make this explana- 
tion unli kely, because Sgr and its tidal tails are —29 kpc 



distant (jSiegel et al.||20lT ). In addition, the Law & Ma- 
jewski ( 2010p model predicts very few high-velocity Sgr 
tidal tail stars in our fields (Fig. H). 

Could these stars be a new substructure in the MW 
halo? We find this explanation unlikely because (1) 
the distribution of high- velocity stars in the dereddened 
CMD (Fig. fTb) is almost identical to the rest of the 
stars in our nelds (which should be dominated by the 
bulge/bar), with rough distances of —5-10 kpc, whereas 
a halo substructure should be more tightly clumped 
around one distance and metallicity; (2) the high- velocity 
stars constitute a large fraction (—10%) of stars in the 
APOGEE fields, which is much larger than would be ex- 
pected for halo substructure; and (3) the high-velocity 
stars are detected in many fields covering a large range 
of longitude, which would imply a truly enormous struc- 
ture. 

We therefore conclude that the high-RV stars are most 
likely members of the bar/bulge. Stars at these veloc- 
ities and in comparable proportions of the total stellar 
population are predicted by the Kazantzidis and MVG 
TV-body models; isolation and examination of these high- 
RV particles in the models reveal their membership in 
the Galactic bar. Furthermore, these model particles are 

f:>uped together on the "far" side (leading edge; see Fig. 
) of the bar, suggesting a coherent spatial/dynamical 
iture, rather than a sampling of stars with coinciden- 
tally similar velocities but on intrinsically unassociated 
orbits. (The lack of high-RV stars in the Besancon model 
may simply be due to their use of highly-smoothed kine- 
matical prescriptions.) 

The nature and significance of the trough at interme- 
diate velocities (Vgsr^ 140-1 80 km s _1 ) is not clear. We 
have isolated the TV-body particles occupying the RV 
trough but could not identify any group property, such 
as large distance, that would preferentially remove the 
corresponding stars from the APOGEE sample. It is 



possible that rather than an intrinsic trough (or drop in 
stellar surface density), the distribution is caused by an 
enhancement in phase space at high-velocity from the 
combination of stars pilin g up at certain phases of their 



combination oi stars pilin g up at certain pnases ol tneir 
orbit (e.g., apogalacticon, Bu reau fc Athanassoula| 1999) 
and projection effects. This explanation is supported 



projection effects. This explanation is supported 
by the structure se e n in the fluid-dynamical model of 
Weiner & Sellwoodl (119991). In their Fi gures 6 and 8, a 
concentration oi high-RV gas exists on the leading edge of 
the bar, with velocities and Galactic positions consistent 
with the high-RV particles "observed" in the Kazantzidis 
and MVG models. This concentration would correspond 
spatially with stars on X2 orbits in the bar, perpendicular 
to the bar's major axis. In this scenario, the observed RV 
trough simply reflects the low count of stars at high RVs 
that are not in this dense dynamical structure. Why this 
property is not reproduced in the TV-body models may be 
due to the resolution of the stellar structures contained 
within the model or the stellar density laws governing 
them. 

The weight of evidence therefore suggests that the best 
explanation for the high-velocity stars is that they are 
members of the MW bar, grouped on particular orbits or 
phases of their orbit. The iV-body models predict that 
high negative velocity stars should be seen at symmetrical 
longitudes in the fourth quadrant. Additional APOGEE 
data, follow-up data in the fourth quadrant, and future, 
more detailed model comparisons will yield insight into 
the nature of these features. 
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